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scribed by the corresponding order cumulant (semi-invariant) of evolution operators 
for the von Neumann equations. For the initial data from the space of sequences of 
trace class operators the existence of a strong and a weak solution of the Cauchy 
problem is proved. We discuss the relationships of this solution both with the s- 
particle statistical operators, which are solutions of the BBGKY hierarchy, and with 
the s-particle correlation operators of quantum systems. 
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1 Introduction 

In the paper we consider the von Neumann hierarchy for correlation operators that de- 
scribes the evolution of quantum correlations of many-particle systems. The necessity 
to describe the evolution of correlations arises in many problems of modern statistical 
mechanics. Among them we refer to such fundamental problems that are challeng- 
ing for mathematics, in particular, the rigorous derivation of quantum kinetic equa- 
tions [2,4-8,15,19,21,27,36], for example, the kinetic equations describing Bose gases 
in the condensate phase [1,16-18,35], and the description of nonequilibrium quantum 
correlations in ultracold Fermi and Bose gases [28] (and references therein). In the pa- 
per we introduce the hierarchy of equations for correlation operators (the von Neumann 
hierarchy) that describes the quantum correlations from the microscopic point of view 
and shows how such dynamics is originated in the dynamics of an infinite-particle sys- 
tem (the BBGKY hierarchy [10,26,33]) and nonequilibrium fluctuations of macroscopic 
characteristics of such systems. 

The aim of the work is to formulate the evolution equations describing correlations 
in quantum many-particle systems with the general type of an interaction potential and 
construct a solution of the corresponding Cauchy problem, then using the constructed 
solution to establish its relationships with the solutions of hierarchies of the evolution 
equations of infinitely many quantum particles. 

We outline the structure of the paper and the main results. In section 2 we introduce 
preliminary facts about the dynamics of quantum systems of non-fixed number of par- 
ticles and deduce the von Neumann hierarchy for correlation operators which gives the 
alternative description of the evolution of states of the many-particle systems. 

In section 3 we define the cumulants (semi-invariants) of evolution operators of the 
von Neumann equation and investigate some of their typical properties. It turned out 
that the concept of cumulants of evolution operators is a basis of the expansions for the 
solutions of various evolution equations of quantum systems of particles, in particular, the 
von Neumann hierarchy for correlation operators and the BBGKY hierarchy for infinite- 
particle systems. 

In section 4 the solution of the initial-value problem for the von Neumann hierarchy 
is constructed and proved that the solution generates a group of nonlinear operators of 
the class C in the space of trace class operators. In this space we state an existence and 
uniqueness theorem of a strong and a weak solutions for such a Cauchy problem. 

In the last section 5 we discuss the relationships of a solution of the von Neumann 
hierarchy for correlation operators both with the s-particle statistical operators, which 
are solutions of the BBGKY hierarchy and with the s-particle correlation operators of 
quantum systems. We also state the general structure of the BBGKY hierarchy generator 
of infinite-particle quantum systems. 
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2 Evolution equations of quantum many- particle sys- 
tems 

2.1 Quantum systems of particles 

We consider a quantum system of a non-fixed (i.e. arbitrary but finite) number of identical 
(spinless) particles with unit mass m — 1 in the space W, v > 1 (in the terminology of 
statistical mechanics it is known as a nonequilibrium grand canonical ensemble [11]). 

oo 

The Hamiltonian of such a system H = H n is a self-adjoint operator with do- 

n=0 

main V{H) = = G f H | i G V{H n ) G n n ^\\ H n^n\\ 2 < oo} C F H , 

n 

oo 

where T-u = 7~L® n is the Fock space over the Hilbert space H (TC° = C). As- 

n=0 

oo 

sume H = L 2 (W) (coordinate representation) then an element ip G Tn = L 2 (W n ) 

n=0 

is a sequence of functions ■?/> = (ijj , V'iOzi), ■ ■ ■ , ipn(<li, ■ ■ ■ , Qn), ■ ■ ■ ) such that H^ll 2 = 
\ipo\ 2 + J2^Li I dqi ■ ■ ■ dq n \ip n (qi, . . . , q n )\ 2 < +oo. On the subspace of infinitely differen- 
tiable functions with compact supports ip n G Ll(W n ) C L 2 (W n ) n-particle Hamiltonian 
H n acts according to the formula (H = 0) 

j_2 n n n 

^n = -y^A^ n + ^ & k \q h ,...,q ik )i; n , (1) 

i=l k=l ii<...<ifc=l 

where <3>( fc ) is a A;-body interaction potential satisfying Kato conditions [30] and h = 2nh 
is a Planck constant. 

An arbitrary observable of the many-particle system A = (A , A±, . . . , A n , . . .) is a 
sequence of self-adjoint operators A n defined on the Fock space Tu- We will consider the 
observables of the system as the elements of the space £(^) of sequences of bounded 
operators with an operator norm [13,33]. 

The continuous linear positive functional on the space of observables is defined by the 
formula 

oo 1 

^-(^o'E^ ( 2 ) 

n=0 

oo 

where Tri v .. in is the partial trace over l,...,n particles, (e°D) = Yl ^f^i,...,nD n is 

71=0 n ' 

a normalizing factor {grand canonical partition function), for which the notation from 
section 5 is used. The sequence D = (I , Di, . . . , D n , . . .) is an infinite sequence of self- 
adjoint positive density operators D n (I is a unit operator) defined on the Fock space Tu- 
This sequence describes the state of a quantum system of non-fixed number of particles. 
The density operators D n , n > 1 (also called statistical operators whose kernels are known 
as density matrices [9]), defined on the n-particle Hilbert space 7i n = 7i® n = L 2 (W /n ), 
we will denote by D n (l, . . . , n). For a system of identical particles described by Maxwell- 
Boltzmann statistics, one has D n (l, . . . , n) = D n (ii, . . . , i n ) if {ii, . . . , i n } G {1, . . . , n}. 

oo 

We consider the states of a system that belong to the space £ 1 (JF W ) = £ 1 (7Y n ) of 

n=0 

sequences / = (/,/!,...,/„,...) of trace class operators f n = / n (l, . . . , n) G £ 1 (7i n ), 
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satisfying the above-mentioned symmetry condition, equipped with the trace norm 

oo oo 

ll/IUl(^K) =53ll/»lltf(M») = Y1 Tr l,..,n|/n(l,---,^)|- 

n=0 n=0 

We will denote by £q the everywhere dense set in £ 1 (jF w ) of finite sequences of degenerate 
operators [30] with infinitely different iable kernels with compact supports. Note that the 
space S^^Ty) contains sequences of operators more general than those determining the 
states of systems of particles. 

For the D G 2} (J-'n) and A G £i(J r n) functional (j2J) exists. One is interpreted as an 
average value (expectation value) of the observable A in the state with density operator 
D (nonequilibrium grand canonical ensemble for the Maxwell-Boltzmann statistics). 

We remark that in the case of a system of a fixed number N of particles (nonequilibrium 
canonical ensemble) the observables and states are the one-component sequences, respec- 
tively, AW = (0, . . . , 0, A N} 0, . . .), £> (A ° = (0, . . . , 0, D N , 0, . . .). Therefore, the formula 
for an average value reduces to 

(AW) = (Tr; v /;.y) ; Tr, N A N D N . 



2.2 The von Neumann equation 

The evolution of all possible states D{t) = (I, Di(t, 1), ... , D n (t, 1, . . . , n), . . .) is described 
by the initial- value problem for a sequence of the von Neumann equations (the quantum 
Liouville equations) [10,31,33] 

jD(t) = -AfD(t), (3) 
D(t)\ t=0 = D(0), (4) 

where for / G 2,\(J-'n) C T>(H) C £ X (JF W ) the von Neumann operator is defined by 

W)n = ~[f n ,H n ] :=- l -(f n H n -HJ n ). (5) 

In the space of sequences of trace-class operators £ 1 (JF- W ) for an abstract initial- value 
problem ([3l)-(jlj) the following theorem is true. 

Theorem 1 ( [24,33]). The solution of initial-value problem (Gj)-^ is determined by the 
formula 

D(t) =U(-t)D(0)U- 1 (-t), (6) 

oo 

where U(—t) = U n (— t) and 

n=0 

U n (-t)=e-^ tH \ 

K\-t) = e* tH », (7) 

Uo(—t) = I is a unit operator. 

For -D(O) G £o(jF w ) c £ l 1 (J r n) it is a strong (classical) solution and for arbitrary 
D(0) G £ 1 (JF W ) it is a weak (generalized) solution. 
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Note that the nature of notations used for unitary groups e ± n tHn is related to 
the correspondence principle between quantum and classical systems (for the later the 
evolution operator for the Liouville equation for the density of the distribution function 
is defined in analogous terms). It is a consequence of the existence of two approaches to 
the description of the evolution of systems based on the description of the evolution in 
framework of observables or states. The evolution operator generated by solution (J6j) for 
/ G 2}{T-h) we will denote by 

U(-t)fU-\-t) :=G(-t)f. (8) 

The following properties of the group Q{—t) follow from the properties of groups (|7j). 

In the space £ 1 (JF W ) the mapping ©:£—>• G(—t)f is an isometric strongly continuous 
group which preserves positivity and self-adjointness of operators. 

For / G £q(JF w ) c T>{— JV) in the sense of the norm convergence S}{Tu) there exists 
a limit [33] by which the infinitesimal generator: —Af = ®^ =0 (—Af n ) of the group of 
evolution operator (jSJ) is determined 

lim±(0(-f)/ - /) = ~{Hf - fH) : = -A/" /, (9) 

where H = (B™ =0 H n is the Hamiltonian (0Q) and the operator: — jfHf — fH) is defined in 
the domain T>(H) C Tn- 

It should be emphasized that the density operator belonging to the space £ X (JF W ) 
describes only finitely many-particle systems , i.e. systems for which the average number 
of particles in the system is finite. Indeed, for the observable of the number of particles 
N = (0, I, 21, . . . , nl, . . .) functional fl2]) has a form 

(N)(t) = (e a ^(0)) ( ; 1 ^-Tr 1) ... )ri+1 L> n+1 (t,l,...,n + l), (10) 

n=0 U - 

and for any sequence D(t) e £ 1 (^-«) we get 

|<A0(*)|< p(0)|U 1( ^)<oo. 

2.3 Derivation of von Neumann hierarchy for correlation oper- 
ators 

Let us represent the state D(t) of a quantum system in the form of cluster expansions 
[10,34] over the new operators g(t) = (0,gi(t, 1), . . . ,g n (t, 1, . . . ,n), . . .) 

D 1 (t,l) = ^ 1 (t J l), 

D 2 {t,l,2)=g 2 {t,l,2)+g 1 {t,l)g 1 {t,2), 



D n (t,Y)= Yl Hg^Xiln^l, (11) 

P:Y=\JXi X,CP 

i 

where the following notations are used: Y = (1, . . . ,n), \Y\ = n is a number of particles 

of the set Y, is a sum over all possible partitions P of the set Y into |P| nonempty 
p 

mutually disjoint subsets X*. 
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It is evidently that, in terms of the sequences of operators g(t), the state of the system 
is described in an equivalent way. The operators g n (t),n > 1 are interpreted as correlation 
operators of a system of particles. 

The evolution of correlation operators is described by the initial- value problem for the 
von Neumann hierarchy for correlation operators 

± 9n (t,Y) = -Af n (Y)g n (t,Y) + (-^ int (Xi,---,X\r\)) Hgm&Xi), (12) 

11 P:Y=(JXi, XiCP 

|P|>1 



g n (t,Y)\ t=(j =g n (0,Y), n>\. (13) 



The von Neumann operator M n {Y) = M n for the system of particles with Hamiltonian 
(CQ) is defined by formula (J5]), 

N M (X 1 ,...,X W ) = E • • E ^lF Zr \z u ...,Z lPl ), (14) 

Z\ C-X"i , -Z|p| C-X"|p| , 
Zi^0 Z| P |^0 

where Yl * s a sum over an subsets Zj C Xj and for k — 1, . . . , n 

^(1,2, ■■•,*) =-![• ,<&«(1,2,...,*)]. (15) 
The simplest examples of von Neumann hierarchy (TT2T) are given by 

j t g 1 {t,l) = -U 1 {l)g 1 {t,l), 

j t g 2 (t, 1,2) = -jV 2 (l, 2)g 2 (t, 1, 2) - .A/g(l, 2) 9l (t, l) 9l (t, 2), 

^ 3 (*, 1, 2, 3) = -jV 3 (l, 2, 3)g 3 (t, 1, 2, 3) + 

+ ( - A/2(l, 2) - ^2(1, 3) - 2, 3))si(t, l)</ 2 (*, 2, 3) + 

+ ( - 2) - ^2(2, 3) - 2, 3))^(t, 2)^(t, 1, 3) + 

+ ( - 3) - A/g(2, 3) - A/g(l, 2, 3))^(t, % 2 (t, 1, 2) + 

-^ ) t (l,2,3)g 1 (t,l)g 1 (t,2)g 1 (t,3). 

We note that, in the case of a two-body interaction potential (k = 2), the von Neumann 
hierarchy (IT21) is simplified. For example, the expression for the operator (73 (t) does not 
contain terms with the operator Af^t- m this case the nonlinear terms in hierarchy (|12j) 
have the form 

E E E (-• A/ ^S^'^))5|x l| (t,^|x 2 |(t,X 2 ). (16) 

P:Y=X 1 \JX 2 heXi «2GX 2 

For classical systems hierarchy (fT2j) with nonlinear terms (fT6|) is an equivalent form of the 
corresponding nonlinear Liouville hierarchy [37] formulated by Green [29]. 
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The von Neumann hierarchy ( TT2l can be formally derived from the sequence of (linear) 
von Neumann equations ([3]) provided that the state of a system is described in terms of 
correlation operators defined by cluster expansions (jTTj) . 

We remark that, in the case of a system of particles satisfying Fermi or Bose statistics, 
cluster expansions (TTTT) and hierarchy (fT2|) have another structure. The analysis of these 
cases will be given in a separate paper. 

3 Cluster expansions of evolution operator of von 
Neumann equation 

3.1 Cluster expansions 

Let us expand the group Q{— t) of evolution operator (JE]) as follows (cluster expansions): 
g n (-t,Y)= £ ]7 ^\ Xi \{t,Xi), n=\Y\>0, (17) 

P:Y=\JXiXiCP 

i 

where * s ^ ne sum over a ^ possible partitions of the set Y = (1, . . . , n) into |P| nonempty 
p 

mutually disjoint subsets X,i C Y. The operators 2t n (t,F) we refer to as the nth-order 
cumulant (semi-invariant) of evolution operators ([8]). 

The simplest examples of cluster expansions (fTT|) have the form 

=2li(t, 1), 

g 2 {-t, 1, 2) = 2l 2 (t, 1, 2) + ^(t, 1) a x (t, 2), 

1, 2, 3) = a 3 (t, 1, 2, 3) + 2l 2 (t, 1, 2)2li(t, 3) + 2l 2 (t, 1, 3) 21^, 2) 
+2l 2 (t, 2, 3) S2ti(t, 1) + 2li(t, l)2li(t, 2)2ti(t, 3). 

Solving previous relations, one obtains 

= &(-*, i), 

2t 2 (t, 1, 2) = 1, 2) - &(-*, l)ft(-t, 2), 

2i 3 (t, i, 2, 3) = g 3 (-t, i, 2, 3) - g x {-t, 3)g 2 (-t, 1, 2) - g^-t, 2)g 2 (-t, 1, 3) - 

-Qi(-t, l)g 2 (-t, 2, 3) + 2!&(-t, l)£i(-t, 2)&(-t, 3). 
In general case the following lemma is true. 
Lemma 1. A solution of recurrence relations [F3\ ) is determined by the expansion 

*n(t,Y)= (-i) |phi (ipi-i)! n *w-*.*o. 

p : y=Ux l x,cp 

i 

n=\Y\> 1, 

where is the sum over all possible partitions of the set Y into |P| nonempty mutually 
p 

disjoint subsets Xi C Y. 
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Proof. Let us consider the linear space of sequences / = (fo, fi(l ),..., f n {l, ■ ■ ■ , n), ... ) 
of operators /„ G £}(H n ) (/ is an operator that multiplies a function by an arbitrary 
number). We introduce in this linear space the tensor *-product [20,34] 

(f*h) lYl (Y) = Y,(f)\z\(Z)h lY \ Zl (Y\Z), (19) 

ZcY 

where h = (0, hi(l), . . . , h n (l, . . . , n), . . .) is a sequence of operators h n G 2}{7i n ) and 

ZcY 

is the sum over all subsets Z of the set Y = (1, . . . , n). 

According to definition (JTSD for the sequence St(t) = (0, 5ti(t, 1), 2t 2 (t, 1, 2), . . . 
. . . , 2ln(t, 1, . . . , n), . . . ) the following equality is true 

E II SW*.**) = (Ex P ,2l(t)) n (F), n = |V| > 1, 

Y>:Y={JXi XiCP 

where Exp^ is defined as the *-exponential mapping, i.e. 

X 

ExpJ = l + E^/^^_/, (20) 

n=1 ' n 

/ = (/ , fx, . . . , f n , . . .) and 1 = (J, 0, 0, . . .) is the unit sequence. 
As a result, we can represent recurrence relations (1171) in the form 

1 + £(-£) =Exp # 2l(t), 

where C/( — t) = (0, i, 1), . . . , {? n (— £, 1, . . . , n), . . .) and the elements of the sequence 
<7(— t) = 1 + Q{— t) = (l,Qi(—t,l),...,Q n (—t,l,...,n),...) are the evolution opera- 
tors ©, (g(-t)f) n (Y) = (U{-t)fU-\-t)) n {Y) = U n {-t,Y)f n {Y)U-\-t,Y) for Y = 
(l,...,n). 

Similarly, defining the mapping Ln* on the sequences h = (0, hi, . . . , h n , . . .) as the 
mapping inverse to Exp„,, i.e. 

00 (—lY 1 - 1 

Ln»(l + /i) =E" h*-y*h , (21) 

n=l n 

one obtains 

E (-l)' Phl (|P|-l)! II Q\x tl (-t,X t )=hn4l + g(-t)) n (Y), 

P:Y=\JXi X,CP 

i 

n=\Y\> 1. 

Hence, relation ([TBI) can be rewritten as 

2l(t) =Ln,(l + (?(-*)), 

and, therefore, expression ( 1T81) is a solution of recurrence relations ( jTTl) . □ 
For systems of classical particles cumulants (|T8l) were introduced in [22]. 
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3.2 Properties of cumulants 

We will now deal with the properties of cumulants (1TB]) . As was proved in Lemma 1 the 
cumulants 2t n (t), n > 1 of evolution operators ([8]) of the von Neumann equations are 
solutions of recurrence relations (jTTj) . i.e. cluster expansions of the group of evolution 
operators ([8]), similar to (TTTj) . 

For the quantum system of non-interacting particles for n > 2 we have: 2l n (t) = 0. 

Indeed, for a non-interacting Maxwell-Boltzmann gas we have: G n (— t, 1, . . . , n) = 

P:Y={JXi XiCPii=l 

-i 

n n 

= 5^(-l)*- 1 B(n J *)(fc-l)!lJ^ 1 (-t l i) = 0. 

fc=i i=i 

where the following equality is used: 

n 

(-l) |Phl (|P| - 1)! = ^(-l^sCn, k)(k - 1)! = 5 nA , (22) 
p : y=u x z k=l 

s(n, k) is the Stirling numbers of the second kind and 5 n \ is a Kroneker symbol. 

In the general case a generator of the nth-order cumulant, n > 2, is an operator 
{—Mint) defined by an n-body interaction potential (|T5|) . According to equality (|22|) for 
the nth-order cumulant, n > 2, in the sense of point-by-point convergence of the space 
^ 1 {'Hn) we have 

\im- t Ql n (t,Y)g n (Y) = £ (-1)I'M(|P| - 1)! ^(-^iC^))^ = 

fe 

= (-^S?(i r ))^(n (23) 

where the operator M^} is given by ( 1T51) . 

For n = 1 the generator of the first-order cumulant is given by 

Hmi^f, Y) - l)g n {Y) = -Af n (Y)g n (Y), 

where Y — (1, . . . , n) and ( — A/" n (Y)) is defined by formula fl5J). 

4 Initial-value problem for the von Neumann hier- 
archy 

4.1 The formula of a solution 

We consider two approaches to the construction of a solution of the von Neumann hi- 
erarchy (fT2l . Since hierarchy (fT2l) has the structure of recurrence equations we deduce 
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that the solution can be constructed by successive integration of the inhomogeneous von 
Neumann equations. Indeed, for solutions of the first two equations we have 

g 1 (t,l) = g 1 (-t,l)g 1 (0,l), 

g 2 (t,l,2)=g 2 (-t,l,2)g 2 (0,l,2) + 

t 

+ J dhg 2 (-t + h,l,2)(- AfSl(l,2))g 1 (-t 1 ,l)g 1 (-t l ,2)g 1 (0,l)g 1 (0,2). (24) 
o 

Then let us consider the second term on the right hand side of (1241) 

t 

J dt x g 2 {-t + t u i, 2) ( - a/2(i, 2))g 1 {-t 1 , 2) = 

o 

t 

/d 
dh— (g 2 {h, i, 2)g 1 (-t 1 , i)&(-ti, 2; 


= g 2 {-t, 1, 2) - g 1 (-t, i)&(-t, 2). (25) 

The operator £ 2 (-t,l,2) - l)C/ a (-t,2) := 2l 2 (t,l,2) in ([25]) is the second-order 

cumulant of evolution operators (jHJ). Formula ( 1251) is an analog of the Duhamel formula, 
which holds rigorously, for example, for bounded interaction potential [3]. 

Making use of the transformations similar to fl25l) . for n > 2 a solution of equations 
(fl2"|) . constructed by the iterations is presented by expressions (|28|) . 

The formula for a solution of the von Neumann hierarchy (fT2l — (fT3"l) can also be (for- 
mally) derived from the solution D(t) = (J, Di(t, 1), . . . , D n (t, 1, . . . , n), . . .) of von Neu- 
mann equations ([HD (Theorem 1) on the bases of cluster expansions (11 II) . 

Indeed, a solution of recurrence equations ffTTj) is defined by the expressions 

g n (t,Y)= (-l) |Phl (|P|-l)' l[D lx ^t,X t ),n>l, (26) 

P:Y=\JXi XiCP 

i 

where is a sum over all possible partitions P of the set Y = (1, . . . , n) into |P| nonempty 
p 

mutually disjoint subsets Xj. If we substitute solution §6§ in expressions (|26j) and taking 
into account cluster expansions (TTTT) for t = 0, we derive 

ff»(t,r)= E (-1) |PK1 (|P|-1)! II E II 9\ Zki \(0,Z ki ). (27) 

P: Y=\JXi X,CP Pi-.X^UZk.Zk.CPi 

i k i 

As a result of the regrouping in expression (1271) the items with similar products of initial 
operators n^cp fl|Xi|(0, X), one obtains 

ffn(t,Y)= E a !P|(*' y p) II ^l(°' X *)' ( 28 ) 

P:Y={JXi X X CP 

i 

where Y = (1, . . . , n), Y P = (Xi, . . . , X| P |) is a set whose elements are |P| subsets Xi C Y 
of partition P : K = [J Xi and ^ is a sum over all possible partitions P of the set Y 

i P:Y=\JX l 
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into |P| nonempty mutually disjoint subsets. Evolution operators 2tip|(£) for every |P| > 1 
in expression ( 128]) are defined by 

2l| P |(t,Yp):= (-l) |P ' hl (iP'|-l)! II S\z k \(-t,Z k ). (29) 

P':Y P =[jZ k Z k CP' 

k 

For |P| > 2 the |P|t/i-order cumulants 2lipi(£) of evolution operators (|S]) of the von Neumann 
equations [25] have similar structure, in contrast to the first-order cumulant. For example, 
for |P| = 1 

2ti(t,lU...Un) =g n (-t,l,...,n), 

for |P| = 2 

2l 2 (Mi U . . . U iiXii^iXxl+i U . . . U i\ Y \) = 

= G\Y\(-t, 1, • • • ,n) - G\x!\{-t,h, ■ ■ • ,«*|x 1 |)^|x 2 |(-^«|x 1 |+i, • • - ,in)> 

where {ii, . . . , iiyi} G {1, . . . , n} and the following notations are used: the symbol (ii U 
. . . U U . . . U Z|y|) denote that the sets . . . , i\x t \} an d {iix^+i, • • • , i\Y\} are 

the connected subsets (clusters, respectively, of \Xi\ and \X 2 \ particles) of a partition of 
the set Y = (1, . . . , |Xl|, |Xl| + 1, . . . , \Y\) into two elements. 
The simplest examples for solution (j28p are given by 

g 1 (t,l) = % 1 (t,l)g 1 (0,l), 

g 2 (t, 1, 2) = 2lx(t, 1 U 2)g 2 {0, 1, 2) + 2l 2 (t, 1, 2)^(0, 1)^(0, 2), 

</ 3 (*, 1, 2, 3) = Sti(t, 1 U 2 U % 3 (0, 1, 2, 3) + 2l 2 (t, 2 U 3, 1)^(0, 1)<? 2 (0, 2, 3) + 

+2l 2 (t, 1 U 3, 2)^(0, 2)0,(0, 1, 3) + 2l 2 (t, 1 U 2, 3)^(0, 3)^(0, 1, 2) + 

+^ 3 (M, 2, 3)^(0, 1)^(0, 2)^(0,3). 

We remark that at the initial instant t — solution fT28l) satisfies initial condition (fT3"l) . 
Indeed, for |P| > 2 according to definitions (jTJ) (W^ bl (0) = 7 is a unit operator) and in 
view of equality fT22|) . we have 

2i |P| (o,yp)= £ (-i)i p 'i- 1 (|p'|-i)!/ = o. 

p' : y P =Uz fe 

k 

4.2 Chaos property 

Let us consider the structure of solution (T28]) for one physically motivated example of the 
initial data; that is to say, if the initial data for Cauchy problem ffT2]) - (fT3]) satisfy the 
chaos property (statistically independent particles) [11], i.e. the sequence of correlation 
operators is the following one-component sequence 

S(O) = (0,^(0,1), 0,0,...). (30) 

In fact, in terms of the sequence -D(O) this condition means that 

D(0) = (l,D 1 (0,l),D 1 (0,l)D a (0,2),...). 
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Making use of relation ( 1261) we obtain initial condition ( 130]) for correlation operators. 
A more general property, namely, decay of correlations for the classical system of hard 
spheres, was considered in [23]. 

For initial data (1301) the formula for solution (|28|) of the initial- value problem (|T2|) - (|T3l) 
is simplified and is reduced to the following formula: 

n 

g n (t,l,...,n) = f& n (t,l,...,n)Y[g 1 (0,i), n>l. (31) 

8=1 

It is clear from ( 1311) that, if at the initial instant there are no correlations in the system, 
the correlations generated by the dynamics of a system are completely governed by the 
cumulants of evolution operators (jHJ). 

In the case of initial data (|30|) solution (l3Tj) of the Cauchy problem for the von Neumann 
hierarchy (|T2|) - (|T3|) may be rewritten in another form. For n = 1, we have 

0i(t,l) = 9li(M)<7i(O,l). 

Then, within the context of the definition of the first-order cumulant, 2li(t), and the 
inverse to it evolution operator 2li(— t), we express the correlation operators g n (t), n>2 
in terms of the one-particle correlation operator gi(t), making use of formula (1311) . 
Finally, for n > 2 formula (l3"Tj) is given by 

n 

g n (t, 1, . . . ,n) = % n (t, 1, . . .,n)Y[gi(t,i), n>2, 

i=l 

where 2l n (t, 1, . . . , n) is the nth-order cumulant of the scattering operators 

n 

Q t (l, ...,n) := Gn(~t, 1, • • • ,n) k), n>l. 

k=l 

The generator of the scattering operator Q t (1, . . . , n) is determined by the operator: 

n n 

-£ E MS (*!,..., u 

k=2 h<...<i k =l 

where N^l ac ^s according to f|T5|) . 

4.3 Properties of a group of nonlinear operators 

On 2} (7i n ) solution (|28|) generates a group of nonlinear operators of the von Neumann 
hierarchy. The properties of this group are described in the following theorem. 

Theorem 2. For g n G £}{H n ), n > 1, the mapping 

t ^ (% t (g)) = J2 2t| P |(t,yp) JI^lTO (32) 

P:Y=\JXi XiCP 
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is a group of nonlinear operators of class Co. In the subspace 2^{7i n ) C 2}{7i n ) the 
infinitesimal generator Af nl (•) of group fifty) is defined by the operator 

{Af nl (g)) n (Y):=-K.(Y)g n (Y)+ £ ( - M**(X U . . . , X m j) J[ g\x A (X,), (33) 

P:Y={JXi, XiCP 
|P|>1 

where the notations are similar to those in [Lfy . 

Proof. Mapping (|32|) is defined for g n e £ 1 (7i n ), n > 1 and the following inequality holds: 

||(2l^))J £W <^ 2 " + V, (34) 
where c := ^ max^ fc|MII £ i (?<|Xi|) - 

i 

Indeed, inasmuch as for g n e £ 1 (?-^ n ) the equality holds [24] 

1h,...,*|0»(-f)</n| = II <7n II 

we have 

ll(a*0/))Jli?(«.)< E E d p 'l-i)! II tailUc«,x 4 ,)< 

|p| |p| 
< E d p lEs(|P|,A;)(A;-l)! < E d p l ^ fcl^ 1 < n !e 2n+1 c n , 
p : y=yXi fe=i p : y=u^i fc=i 

where s(|P|,A;) are the Stirling numbers of the second kind. That is, (%k(g)) £ ^(Wn) 
for arbitrary iGl 1 and n > 1. 

We can now formulate the group property of the one-parametric family of nonlinear 
operators which are defined by (1321) . i.e. 

9t tl (2l* 3 (^)=Stt a (2t* 1 (p))=St tl+ta U). 

Indeed, for g n G -C^W^), n > 1 and for any ti, t 2 G R 1 , according to ( 1281) and ( lT8i) . we 
have 

(%0Ms)))n( y ) = E %i(*i^p)II E %i(*2,{^} Pi ) J] = 



e e (-i) |p ' hi (ip'i-i)! n ^iHi.q 

P:y=U^iP':Y P =UQj QjCP' 



x n e e (-i) ip ' m (ip:i-i)< n w-^) n 

^cpp 1 :jr 1 =uz I . p ; :{Xl}p . =Uflfc . ^cp; z^cPi 



where {Xj} Pi = (Zi, . . . , Zpj) is a set whose elements are |Pj| subsets C Xj of the 

partition Pj : — \}Z[.. Having collected the items at identical products of the initial 

k 
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data <7n(0), n > 1, and taking into account the group property of the evolution operators 
U^(-t), n > 1 (CD, we obtain 

0M^)))„On= E E (-l) |P ' hl (|P'|-l)' II G\Q k \(-ti-t2,Q k )x 

P:y=UX lP ' : y p=UQfc QfcCP ' 

» k 

x n »w = e ^\(ti+t 2 ,Y P ) n ^pQ) = (2t tl+t2 (^) n (y). 

XiCP P:Y=\JXi XiCP 

i 

Similarly, we establish 

The strong continuity property of the group 2l t (g(0)) over the parameter t G IR 1 is a 
consequence of the strong continuity of group ([8]) of the von Neumann equations [13]. 
Indeed, according to identity (1221) the following equality holds 

E E (- 1 )' P ' hl (l p 'l - !)' II W*) = 9n(Y). 

P:Y=UX lp i :Yp =\JZ k ^CP 

i k 

Therefore, for g n G £o(7i„) c ^C^n); ^ > 1, we have 



fell £ £ (-l) |P ' hl (|P'|-l)! II G\z k \(-t,Z k ) \{g\x^)- gn (Y) 

P.Y=\JX,p':Y P =(JZ k Z k CP' ^iCP 



< e e (i p 'i-i)!jim|| n n^iw-n^i(^ 

p : y=U^p':y P =u^ ^ z k cp' **cp ^cp 



In view of the the fact that group £/ n (— (jBD (Theorem 1) is a strong continuous, i.e. in 
the sense of the norm convergence £}{7i n ) there exists the limit 

\\m(g n (-t)g n - g n ) = 0, 

which implies that, for mutually disjoint subsets X$ C Y, the following equality is also 
valid: 

lim( Yl G\z k \(~t, Zk)9n - 9n) = 0. 

Z k CP' 

For g n G £j(W n ) C £ 1 (7i n ), we finally have 

\im\\(%(g)) n - g n \\ £ i {Hn) = 0. 

We will now construct the generator A/" ni (-) of group (131Z1) . Taking into account that for 
g n G ^(Tin) C P(A/"™'(-)) equality holds, let us differentiate the group (2(i(<?)) ip n for 
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all ip n G V(H n ) C Tin- According to equality ( |23i) for the |P|i/i-order cumulant, |P| > 2, 
we obtain 

gmis2t,p,(t,y P )(7 n V» = E (- 1 )' P ' M (I P 'I " X ) ! E (--% fc |(^))<?nV>n = 



p':Y P =UZ fe 



ZfcCP' 



.( £ \Zr\) 



= E • ■ E ("^r 1 ' (^l,...,^|P|))^n= (-AT m *(r P ))^n, 

ZidXi, Z\ P \(ZX\p\, 

where Yp = (Xl, . . . , X\ P \) is a set whose elements are |P| subsets X{ C Y of partition 
P : K = (J Xj and the operator Af^t is given by formula (TT5T) . 

i 

Therefore, for group fl32l) we derive 



! im K( 2l *^))n-^)^ = !™7( E 2i|pi(^^) n^i(^)-^( y ))^ 

i 

= \im j(Vl 1 (t,Y)g n -g n )ij n + E j™ T 2 ^' F p) II = 



P:Y=UXi 
|P|>1 



X,CP 



(-J^„)(y)i+ E (-^ int (^p)) n^iW^- ( 35 ) 



p : y=U^i, 
[p|>i 



XCP 



Then in view of equality fl35l) and the proof of the theorem 2 for g n G £Q(7i n ) C 
T>{N™ {•)) C £ 1 (7^ n ), n > 1, in the sense of the norm convergence in £ 1 (?-^ n ), we finally 
have 



lim 

t->o 



■((2l t (,)) n -,„)-(^(,)) ? 



0. 



where M nl {-) is given by formula ([33 



□ 



4.4 The uniqueness and existence theorem 

For abstract initial- value problem (jT2l) - (jT3l in the space £}{7i n ) of trace class operators 
the following theorem holds. 

Theorem 3. The solution of initial-value problem [W) - U~3) for the von Neumann hier- 
archy ^^) is determined by formula For g n (0) G £j(7Y n ) C 2}{7i n ) it is a strong 
(classical) solution and for arbitrary initial data g n {0) G £ 1 (7i n ) it is a weak (generalized) 
solution. 

Proof. According to theorem 2 for initial data g n (0) G £q(7Y„) C £ 1 (?-^ n ), n > 1, sequence 
(1281) is a strong solution of initial- value problem (fT2l - (jT3l) . 
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Let us show that in the general case g(0) G expansions ( l28i) give a weak solution 

of the initial- value problem for the von Neumann hierarchy ( jT2i) . Consider the functional 

(^n,9n(t)) := Tti,..., n Pnflw(t), (36) 

where G £o(W„) are degenerate bounded operators with infinitely times different iable 
kernels with compact supports. The operator g n (t) is defined by (1281) for the arbitrary 
initial data gk(0) G H 1 ^*;), = 1, . . . , n. According to estimate (I34I) for g n G and 
f n G £ ( / ^n) functional (155]) exists. 

We transform functional (l3T)j) as follows 



(Vn,^n(*)) = E (fn,%\p\(t,Y P ) Yl 9\Xi\(0,Xi 



P:Y=(JX t XiCP 

i 

= e e (-i) |p ' hi dp'i-i)!( n w^k, n^i(o ; ^))>(37) 

p : y=UX lP ' : y p= yz fc Z fc CP' ^iCP 

» fc 

where the group of operators Q n {t) is adjoint to the group Q n {—t) in the sense of functional 
flU. 

For g n (0) G £ 1 (H n ) and <p n G £o(7i n ) within the context of the theorem 2 we have 

lfm ({q{Gn{t)Vn ~ V 3 n),5'n(0)) - (MnVn, 5^(0))) =0, 

and, therefore, it holds that 

^ Z fe CP' ^CP 

II 0i**iH; z *) II 

Z iC P' Z fc CP' *iCP 



Then, for representation (1371) of functional (|36|) . one obtains 

J t (<Pn,gn(t)) = (K<Pn,9n(t))+ E E ( _1 ) ( l P ' I ~ 1)' 

P:Y=U^,p':y P =UX fc 

1 fe 
|P|>1 

x( ^2 Mz^z^n, n f?iz fc i(-t;^) n 9&*\{o,Xi 

ZiCP' Z k CP' x ^ p 

|P|>1 

where the operator J\f mt (Xx, . . . , X\ P \) is defined by (TTi|) and (TT5|) . 
For functional (1361) we finally have 

= ((K<Pn)(Y),g n (t,Y))+ (j^ nt (X 1 ,...,X ]P \) l p n , J] (38) 

p : y=U^i, ^cp 
|p]>i 
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Equality ( 1381) means that for arbitrary initial data g n (0) G ^(TCn), n > 1, the weak 
solution of the initial-value problem of the von Neumann hierarchy (|T2l)-(jT3l) is determined 
by formula ([28]). □ 



5 BBGKY hierarchy 



5.1 Nonequilibrium grand canonical ensemble 

As we have seen above the two equivalent approaches to the the description of the state 
evolution of quantum many-particle systems were formulated, namely, both on the basis 
of von Neumann equations ([3]) for the statistical operators D{t) and of the von Neumann 
hierarchy ffT2]) for the correlation operators g(t). For the system of a finite average number 
of particles there exists another possibility to describe the evolution of states, namely, by 
sequences of s-particle statistical operators that satisfy the BBGKY hierarchy [10]. 

Traditionally such a hierarchy is deduced on the basis of solutions of the von Neumann 
equations (the nonequilibrium grand canonical ensemble [11,25,33] or the canonical en- 
semble [2, 10, 15]) in the space of sequences of trace class operators. 

The sequence F(t) = (I, Fi(t, 1), . . . , F s (t, 1, . . . , s), . . .) of s-particle statistical opera- 
tors F 8 (t, 1, . . . , s), s > 1 can be defined in the framework of the sequence of operators 
D{t) = (I, Di(t, 1), . . . , D n (t, 1, . . . , n), . . .) (the operator D n {t) being regarded as density 
operator ([6]) of the n-particle system) by expressions [33] 

F s (t, 1, . . . , s) = {e a D(0)) Q 1 -jTr s+li ... jS+n L> s+n (t, l,...,s + n), s>l, (39) 

n=0 

oo 

where (e a Z)(0)) = 1 + ^T^i,...,nD n (0, 1, . . . ,n) is a partition function (see definition 

n=l 

of functional ([2])). For D(0) G ^(JF-^) series ( 1391 converges. 

If we describe the states of a quantum system of particles in the framework of correla- 
tion operators g(t) the s-particle statistical operators, that are solutions of the BBGKY 
hierarchy, are defined by the expansion 

oo ^ 

F s (t, 1, . . . , s) = — Tr s+lt _ >s+n g 1+n (t, 1 U . . . U s, s + 1, . . . , s + n), s > 1, (40) 

n=0 n ' 

where gi +n (t, lU...Us,s + l,...,s + n), n > 0, are correlation operators (|28|) that 
satisfy the von Neumann hierarchy (|T2|) for a system consisting both from particles and 
the cluster of s particles (1 U . . . U s is a notation as to formula (1291 ). 

Expansion (1401) can be derived from (I3"9"|) as a result of the following representation for 
right hand side of expansion fl39|) 



F s (t,l,...,s) 



oo 



S + l,...,S+1l 



E 



-i)i-'-'(iPi-i)! n^,i(«-f.),f>i, 



n=0 



P:{lU...Us,s+l,...,s+n}=|J^ 



(41) 
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where is a sum over all possible partitions P of the set {lU...Us,s + l,...,s + n} 
p 

into |P| nonempty mutually disjoint subsets Xj. If D(t) G 2}(Tu) series ( 14"T1) converges. 

We remark that expansion (1391) can be defined for more general class of operators then 
from the space 

To prove representation f|4T]) we will first introduce some necessary facts. We will use 
the notations of lemma 1. According to the definition of ^-product (1T91) for the sequence 
/ = (/o, /2(1, 2), • • • , / n (l, • • • , n), . . . ) of elements /„ G 2}(H n ) the mapping Exp„ 
is defined by series ( 1201) . The mapping Ln* that inverse to Exp^ is defined by series (12"T1) . 

For / = (/„, /i, ...,/„,. • •), /„ G 2}{U n ) we define the mapping 9 X : / -> Di/ by 

(f i/)n(l, . . . , n) := /„+i(l, . . . , n, n + 1), n > 0, 

and for arbitrary set K = (1, . . . , s) we define the linear mapping 0i . . . D a : f — > Vi . . . D s f 
by 

(?)i...d a /) n := / s+n (l,...,s + n). (42) 

We note that for sequences f l = (f£, f{, . . . , ft, . . .), ft G 2}(H n ) and f 2 = 
Ul fl ...,ft,.. ft G £\H n ) the following identity holds [34] 

0i(/ 1 */ 2 ) = c) 1 / 1 */ 2 + / 1 *c) 1 / 2 . 

Further, since Yp = (Xi, . . . , X\ P \) then Y\ — (1 U . . . U s) (cluster of s particles) is one 
element (|YjJ = 1) of the partition P (|P| = 1), we introduce the mapping X) Yl '■ f — > fyj/, 
as follows 

(D Y J)n(l,---,n) = f n+1 (Y 1 ,l,...,n,), n>0. (43) 

Then for arbitrary / = (/ , /i, ...,/„,. . .), /„ G £}(H n ) according to definition (12171) of 
the mapping Exp^ the following equality holds 

c) 1 (ExpJ) = c) 1 /*ExpJ. (44) 

For f n G £ 1 (H n ) an analog of the annihilation operator is defined 

(a/),(l, ...,s):= Tr s+1 f s+1 (l, ...,s,s + l), (45) 

and, therefore 

oo ^ 

(e a f) s (l, ...,s) = y~] — Tr s+1 s+n / s+n (l, . . . , s, s + 1, . . . , s + n), s>0. 

n! 

n=0 

Using previous definitions and ( ITOj) for sequences ft = (ft, ft, . . . , ft, . . .),ft G 2}{7i n ) 
and ft = (ft, ft, ...,ft,.. .), ft G Z l (H n ) we have [34] 

(e a (ft * / 2 ))o = (e7 1 )o(e7 2 )o- (46) 

From equality (Jl6|) we deduce that expressions (139|) and (jlQ]) can, respectively, be rewritten 

as 

F s (t,Y) = (e a D(0)) 1 (e a D Yl D(t)) 

and 

F.(t,y)= (e a fy^W) , 

where K = (1, . . . , s). 

Hence, in view of equalities (1TT1) . (12T1) to derive expressions (H0|) we have to prove the 
following lemma. 
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Lemma 2. For f = (/o, f\, . . . , f n , . . .) and f n G £}(7i n ) the following identity holds 

(e Q ExpJ) - 1 (e a c)y 1 (ExpJ)) = (e n c)yJ)o. (47) 
Proof. Indeed, using equalities (j4"51) . (jHJ) and (JlSj) . we obtain 

(e a ExpJ) 1 (e a 0y 1 (ExpJ)) = (e n Exp J) 1 (e"(Exp J * a n /)) = 
= (e a ExpJ)o 1 (e a ExpJ)o(e°c)y 1 /) = (e^yj) . 

□ 

5.2 On a solution of the BBGKY hierarchy 

In this subsection we will turn to the solution of the initial-value problem of the BBGKY 

oo 

hierarchy. In the space £^ (JF W ) = a n £ 1 (?i n ), where a > 1 is a real number, we consider 

n=0 

the following initial-value problem for the BBGKY hierarchy for quantum systems of 
particles obeying Maxwell-Boltzmann statistics 

jF s (t) = -KF s (t) + 

oo 

+ E ^T Tr >-' E ( " ^nt l+n) (Z, 8 + 1,..., a + n))F s+n (t), (48) 

n=l ' ZcY, 

F s (t) | t=0 = F,(0), s > 1, (49) 

where Y = (1, . . . , s) and the operator M^} is defined on £„ (J^) C £} a {T-H) by formula 
( |T5l) . The equation from hierarchy (l4~8i) for s = 1 has the following transparent form: 

-F^t) = -M x Fi(t) + ^ T ^,.,n+i ( - <i; +1) (l, 2, . . . , n + l))F n+1 (t). 

n=l 

In the space ^(JF^) there is an equivalent representation for the generator of the BBGKY 
hierarchy: 

oo 1 

e«(-A/>-° = -N + J] -[. . . [AT, a], . . . , a] 

71=1 • 

n— times 

that follows from fl39l) . For a two-body interaction potential one is reduced to the form [24]: 
-M + [Af, a]. 

We remark that in terms of the s-particle density matrix (marginal distribution) 
F s (t, qi, . . . , q s ; q v . . . , q s ) that are kernels of the s-particle density operators F s (t), for 
a two-body interaction potential (see ([T]) for k = 2) the evolution operator (T48"l) takes the 
canonical form of the quantum BBGKY hierarchy [10] 

2 



ih dt Fs< " t] Qh ■ ■ ■ ' qs] q ' 1 ' ■ ■ ■ ' ^ = ( ~ y E( A ^ - A <?;) + 



i=l 

s 

+ E ( $ (* - - - 9i, • • • > q'i, ■ ■ ■ , q's) + 



i<j=i 

s 



+ E / dqs+l (®( qi ~ qs+1 "> ~ _ Qs+i))F s+ i{t; q u ... q s , q s+1 ; q[,...,q s , q s+x ) 
i=i J 
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For the solution of the BBGKY hierarchy the following theorem is true [25]. 

Theorem 4. If F(0) G and a > e, then for t G 1R 1 there exists a unique solution 

of initial-value problem ffify-ffllfy given by 

F s (t,l,...,s) = 

= ^ — Tr s+1> ... )S+n 2li +n (t, Y"i, s + 1, . . . , s + n)F s+n (0, l,...,s + n), s > 1, (50) 

n=0 n ' 

where for n > 

a 1+n (t,y 1>s + i,..., s + n)= ^ (— i) |:p| " :L (|P| — i)! n 

P:{Y u x\Y}=\jXi x iC p 

zs i/ie (1 + n)th-order cumulant of operators X \ K = {s + 1, . . . , s + n} and Y\ = 
{lU...Us}. 

For initial data F(0) e £„ C S^J-h) a strong solution and for arbitrary initial 
data from the space £^,(^ 7 ^) it is a weak solution. 

The condition a > e guarantees the convergence of series (1501 and means that the 
average number of particles ( JTUl) is finite: (N) < oc/e. This fact follows if to renormalize 
a sequence (150]) in such a way: -F s (t) = (N) s F s (t). 

For arbitrary F(0) G the average number of particles (expectation value (jTOl) 

expressed in terms of s-particle operators) 

(N)(t)=Tr 1 F 1 (t,l) (51) 

in state (!50l) is finite, in fact, 

|(iV)(t)|<c Q ||F(0)|| Ww) <oo, 

where c Q = e 2 (l — -) _1 is a constant. We emphasize the difference between finite and 
infinite systems with non-fixed number of particles. To describe an infinite particle system 
we have to construct a solution of initial-value problem (I48I) - (I49I) in more general spaces 
than £^,(JF W ), for example, in the space of sequences of bounded operators to which the 
equilibrium states belong [20,34]. 

We remark that the formula for solution fl50|) can be directly derived from solution 
( |28l) of von Neumann hierarchy ( fT2l [26]. In papers [2,27,32,35] a solution of initial- 
value problem ( l48l) -(T49l) is represented as the perturbation (iteration) series, which for a 
two-body interaction potential has the form 

t t n -i 

OO „ „ s 

F s(t) = J2 M*i- - / <Tr s+w+n &(-mi)^ (-A/™(i 1 ,s+i))g 8+1 (-t 1 +t 2 )x 

n=0 oo il=1 

s+n— 1 

n ) E (-^nt(^s + n))g s+n (-t n )F s+n (0). (52) 
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In the space expansion ( )48i) is equivalent to this iteration series. This follows 

from the validity of analogs of the Duhamel formula for cumulants of evolution operators 
of the von Neumann equations. For the second-order cumulant an analog of the Duhamel 
formula has the form (125]) and in the general case the following formula takes place 

2li +n (t, Yi, s + 1, . . . , s + n) = 

t tn-l 

= Jdt x ... J dt n ]J g^-t + txM) E E E (-JV2(*i>™0) x 

q o kjeli ii<ji£li heii mieji 

xn^lK + ^y E E E (-^{h,m 2 )) x... 

fc2 ££2 *2 <J2 £2 &2 rn-z GJ2 

•••X ] [ £|fc„|(-tn-l+t n , fcn) E E E i~^int ( l n, m n ))Q s+n (-t n , 1, . . . , S+Tl), 

where 2i = {Yi, s + 1, . . . , s + n}, Z n = {?„_i U j n _i} U J n _i \ {i n -i,i n -i}- 

We remark that for classical systems of particles the first few terms of the cumulant 

expansion (|50|) were considered in [12]. 

In [25] we discuss other possible representations of a solution of the BBGKY hierarchy 

in the space 2} a {Tu)- 

5.3 Correlation operators of infinite-particle systems 

Correlation operators f[2"gj) may be employed to directly calculate the macroscopic values 
of a system, in particular, fluctuations characterized by the average values of the square 
deviations of observables from its average values. For example, for an additive-type ob- 

n 

servable a = (a , ai(l), • • • , X) a i(^)> • • •) from the formula for expectation value (j2J) we 

i=i 

derive the formula for fluctuations (the dispersion of an additive-type observable) [10] 

((a-(a)(t)) 2 )(t) = 

= Tr x (a?(l) - (a) 2 (t))F!(t, 1) + Tr 1 , 2 a 1 (l)a 1 (2) (F 2 (t, 1, 2) - F 1 (t, l)iq(t, 2)). 

Therefore, the dispersion of the additive-type observable is defined not directly through 
the solutions of the BBGKY hierarchy but by the following correlation operators: 
F 2 (t, 1,2) — Fi(t, l)Fi(t, 2) = G2(t, 1,2) or in the general case by the s-particle corre- 
lation operators 

G.{t,l,..,s):= E (-1)' PM (|P|-1)! II^I^)- ( 53 ) 

P:{l,...,s}=yXi *iCP 

i 

The s-particle correlation operators G s (t), s > 1 can be expressed in terms of correlation 
operators fl28|) by the formula 

G s (t, l,...,s) = E~| ^+1 s +nff>+n(M.---,s + n), s > 1, (54) 

n=0 

where g s+n (t, 1, . . . , s + n) is a solution of the von Neumann hierarchy (fT2l) . 
To derive expression (1541) we state the following lemma. 
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Lemma 3. Let f = (/ , fi, ■ ■ ■ , f n ) and f n G £ 1 ('H n ), then the equality holds 

(e a ExpJ)^e a ExpJ = Exp*e a f. 
Proof. Indeed, using equality fj46|) and the equality 

Oi . . . D s ExpJ = Exp J * ^ D x . . .5|Xi|/ * • • • *^|X|p|_i|+i • • -Q\X\t\\f, 

that follows from fj42|) . (j44l) . one obtains (here y = (1, . . . , s)) 

(e a ExpJ) 1 (e a ExpJ)| y |(F) = (e°Exp J)o\e% . . . D s (Exp J)) = 

= (e a Exp J) 1 (e a (ExpJ * ]T 0, . . . V; / *••• *\x w _ L \+i ■ ■ ■ * |x |P| |/) 

= (e a ExpJ) 1 (e a Ex P J) ('e a ^ a . . . X>\ Xl \f *■■■* ^x^+i ■ ■ -5|x |P| |/ 

^ p : y=UX 

= E n( e ^---W)o= (Exp,eV) |y| (y). 
p : y=UX ^cp 

i 

□ 

We now derive representation ( |5~^1) for s-particle correlation operators. Using equality 
(ITT]) in terms of the mapping Exp„, ( l20l) . representation ( 1391) for s-particle statistical 
operators can be rewritten in the form 

F(t) = (e°Exp^(0)) 1 e a Exp^(t). 

In terms of the mapping Exp* formula ( 1531) has the form 

F(t)=Ex^G{t). 

Further, according to previous formula and Lemma 3 we have 

Exp,G(t) = (e a Exp^) VExp^(t) = Exp*e a g{t), 

and, therefore, we finally derive 

G(t) = e a g(t) (55) 

or in component-wise form (j54p . 

For chaos initial data (1301) obeying Maxwell-Boltzmann statistics according to (1551) 
we have Gi(0) = gi(0). Using solution (131]) of the von Neumann hierarchy (|T2l) the 
expansion for a solution of the initial-value problem for s-particle correlation operators 
can be represented as follows 

x j s+n 

G s (t, l,...,s)=) j — Tr s+ i i ... )S+n 2t s+n (t, 1, . . . , s + n) JJ Gt(0,i), s > 1. 

n=0 ' i=l 

It should be noted that sequence (l55j) is a solution of the nonlinear BBGKY hierarchy 
for s-particle correlation operators [29] which describes the dynamics of correlations of 
infinite-particle systems. 
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